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ABSTRACT 

Several attempts have been made in the past to assess the expected number of 
exoplanetary transits that the Gaia space mission will detect. In this Letter we 
use the updated design of Gaia and its expected performance, and apply recent 
empirical statistical procedures to provide a new assessment. Depending on the 
extent of the follow-up effort that will be devoted, we expect Gaia to detect a 
few hundreds to a few thousands transiting exoplanets. 

Subject headings: Methods: statistical - Planets and satellites: detection - Sur- 
veys - Techniques: photometric 



1. Introduction 



Gaia is a planned European Space Agency (ESA) mission, scheduled to be launched 
at 2013. It will perform an all-sky astro metric and spectrophotometric survey of point-like 
objects between 6th and 20th magnitude. The primary goal of the telescope is to explore the 
formation, dynamical, chemical and star-formation evolution of the Milky Way galaxy. The 
main science product of Gaia will be high precision astrometry, backed with photometry 
and spectroscopy. It will observe about 1 billion stars, a few mill ion galaxies, half a million 
quasars, and a few hundred thousands asteroids (ILindegrenl 120 1 Ol ) . 

Gaia will operate in a Lissajous-type orbit, around the L2 point of the Sun-Earth 
system, about 1.5 million kilometers from Earth in the anti-Sun direction. It will have 
a dual telescope, with a common structure and common focal plane. During its five-year 
operational lifetime, the spacecraft will continuously spin around its axis, with a constant 
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speed of 60arcsecs _1 . As a result, during a period of six hours, the two astrometric fields 
of view will scan all objects located along the great circle perpendicular to the spin axis. 
As a result of the basic angle of 106.5° separating the astrometric fields of view on the sky, 
objects will transit the two fields of view with a delay of 106.5 minutes. Due to the spin 
motion of six-hour period, and a 63-day-period precession, the scanning law will be peculiar 
and irregular. This scanning law will result in a total of 70 measurements on average for 
each celestial object Gaia will observe ( Ide Bruijndl2012l ). 



Gaia will provide photometry in several passbands, the widest of which will be a 'white' 
passband dubbed G, centered on Ao = 673 nm, with a width of AA = 440 nm. In what follows 
we use the apparent G magnitude as approximately equal to the apparent V magnitude. 
One can expect a milli-magnitude (mmag) precision in the G band for most of the objects 
Gaia will observe, down to 14th-16t h G magnitude, and 10 mmag at the worst case of 19th 
magnitude objects (IJordi et al.ll2010l ). The exact limiting ma gnitude for a lm mag precision 
depends on instrumental factors which are not yet "frozen" (jde Bruijndl2012l ). 



The 1 mmag precision of Gaia photometry naturally raises the question of whether it 
can be used to detect exoplanetary transits. While 1 mmag precision is nominally more 
than sufficient for the detection of Jovian transiting planets, the low cadence and the small 
number of measurements make the feasibility of this detection a non-trivial question. In 
the literature, there have been several conflicting estimates as to the number of transits 
detectable by Gaia, based on different assumptions. 



Hpgi (120021 ) estimated the expected number of transit detec tions by G aia at 250, for 



long-period planets, with orbital radius of 0.5 AU. According to iHogl (120021 ). the expected 
yield of hot Jupiters (HJs) and very hot J upiter s (VHJs), with orbital radii of 0.1 — 0.01 AU 
was 6250 and 625000 respectively iHogl ( 120021 ) based these estimates on some general as- 
sumptions. First, the p lanets frequency was ap proximated as ~ 1%, w hich a lready proved as 



overestimated for HJs faeattv fc Gaudill2008l 



, hereafter BG). Second, iHogi did not account 



for the stellar density and its variation due to galactic structure, and neglected extinction 
by dust. He assumed that the number of stars Gaia will observe with 1 mmag precision (up 
to magnitude G = 15.5) is 10 9 . Finally, he assumed that a transit detection can be made 
with only one transit observation per system, counting on Gaia astrometric observations to 
complement the transit observation. 



Robichonl (120021 ) performed transit simulations with the assumed Gaia photometry to 



estimate the number of detections, and concluded that Gaia will detect between 4000 and 
40000 transiting Jupiter-like planets. He assumed that the number of individual observations 
per sta r would be between 1 00 and 300 with an average of 130. He also used a galactic 
model (IHaywood et ajj 119971 ). and derived the probability distribution of the number of 
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observations during transits using Gaia scanning law. We suspect that these predictions are 
overestimated, due to the fact that the currently planned scanning law of G aia implies an 



average of 70 measurements per star, and not 130. We will refer to iRobichonl 's estimates in 
more detail in section El 

In this short Letter we present a new estimate, one that we believe is more realistic than 
the previous ones. It is based on the methodology of BG, which includes broad assumptions 
about the galactic structure, as well as implicit assumptions about the geometric transit 
probability and the effects of stellar variability, that rely on statistics from completed transit 
surveys. The field of planetary transits of HJs and VHJs seems to have come to a certain 
maturity from which we can draw some statistical assumptions. We feel it is not yet the 
case for transits of smaller planets ('Saturnian' and 'Neptunian'), and we therefore do not 
address these issues here. 

Providing more accurate and up-to-date predictions of the number of transiting plan- 
ets detectable by Gaia is important for the ongoing effort of developing the Gaia analysis 
pipeline. Moreover, detection of planetary transits requires considerable follow-up effort, 
which is also a reason for having a reliable estimate of this number. We took upon ourselves 
to provide a somewhat more rigorous analysis, based on empirical statistics, simply because 
the estimates in the literature are too varied and inconsistent. At the time they were made, 
the field of transit surveys was too young, and much had yet to be learned. The time has 
come to provide a more decisive estimate, based on the current more evolved understanding 
of the problem. 



2. Predicting the transit yield 

BG presented a statistical methodology to predict the yield of transiting planets from 
photometric surveys. Their method takes into consideration the frequency of short period 
planets, variations in the stellar density due to the galactic structure and also corrects for 
the extinction by dust. 

Following the procedure suggested by BG, the average number of exoplanets that a 
photometric survey can detect is estimated by the product of the probability to detect a 
transiting planet, the frequency of transiting planets and the local stellar mass function. 
Obviously, this product should also be integrated over mass, distance and the field of view: 

-p*(r,l,b)r cosb——— — P det (M,r, R p ,p), (1) 



dR p dp dM dr dl db dM dR p dp 
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p*(r, /, b) is the local stellar density as a function of heliocentric galactic coordinates (r, /, b), 
dn/dM is the present day mass function in the solar neighborhood, and d is the 

frequency of transiting planets (the probability that a given star harbors a transiting planet 
with radius R v and period p). Pdct is the transit detection probability, assuming there is 
indeed a transiting planet around the examined star. 

Following BG we consider detection of transits of Very Hot Jupiters (VH Js) with orbital 
periods of 1 — 3 days, and Hot Jupiters (H Js) with periods of 3 — 5 days. For the probability 
that a star will harbor a transiting planet with radius R p and orbital period p we assume 
the same form that BG used: 



d 2 f(R p ,p) 
dR p dp 



k(p)f(p)5(R p -R' 



(2) 



The normalization factor k(p) is an empirical number which can be deduced from com- 
pleted surveys. It encapsulates many factors which are common to transit surveys looking 
for the same range of periods and more or less the s ame stellar populat ions. BG proposed 
to use a normalization factor based on the results of iGould et al.l (120061 ). who analyzed the 
statistics of the OGLE transit surveys. They suggested a value of k(p) = 1/690 for VHJs, 
and k(p) = 1/310 for HJs, and a locally uniform distribution of the period in the specified 
interval, f(p) = 1. 



This planet frequency ( Gould et al. 20061) is compatible with the frequency which is 
implied by Kepler results ([Howard et al.ll201ll ). Furthermore, Gaza's low cadence makes it 
more similar to ground-based surveys, rather than to high-cadence space surveys. Moreover, 
CoRoT and Kepler results are not yet complete, since they are still operating. 

To account for the stellar densi t y in t he solar neighborhood, we used the Present Day 
Mass Function (PDMF). lReid et al.l (120021 ) used data from Palomar /Michigan State Univer- 
sity survey, together with the Hipparcos dataset, to derive the PDMF: 



dn 
dM 



fcnormd^)- 1 ' 35 for 0. 1 < M/M < 1 



h ( M ^-5.2 
'"norm V M& ) 



for 1 < M/M c 



(3) 







where again we adopted the normalization suggested by BG, k xi 



0.02124 pc" 



In order to con yert absolut e magn itudes to masses, we used the mass-luminosity rela- 
tions, as derived by iReid et al.l ( 120021 ). For lower main-sequence stars we used the mass- 
luminosity relation: 



logM = 10~ 3 x (0.3 + 1.87M V + 7.6141/2 



1.698M 3 



0.06096M£), 



(4) 
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and for the upper main- sequence stars, 



log M = 0.477 - 0.135M y + 1.228 x 10~ 2 Mf> - 6.734 x 10" 4 M, 



Vi 



(5) 



where My is the absolute visual magnitude of the star. The boundary between the calibra- 
tions is set at My = 10. 

The last step of the procedure is integration of the stellar density over the entire field 
of view: 



N, 



dot 



p* (r, I , b) r 2 cos b db dl dr. 



(6) 



We now had to account for the variation in the stellar density due to the galactic 
structure, and also the effect of e xtinction due to int e rstella r dust. We incorporated into our 
calculations the galactic model of iBahcall fc Soneiral (119801 ). which is a simplified model that 
depends only on the distance from the galactic plane (z) and the disk-projected distance 
from the galactic center (d): 

d \z\ 



h, 



d.* 



2.5 kpc is the scale length of the disk. The scale height, H(My) 



(7) 

depends on the 



absolute magnitude, and BG use the dependence: 

90 for My < 2 

H(M V ) = { 90 + 200 ( M ^= 2 ) for 2 < My < 5 



290 



for M v > 5 



To account for the interstellar dust extinction, we used the expression suggested by 
Bahcall fc Soneiral (1l980l ) for obscuration by dust. In their model the obscuration depends 
on the heliocentric distance and the galactic latitude: 



A{r) 





0.17 (1.2 + tan |6|) [1 - exp( 



\b\ > 50° 
-^)] esc |6| \b\ <50° ' 



(9) 



where h is a typical scale height with respect to the galactic plane, that is approximated by 
lOOpc feahcall fc Soneiral Il980h . 



The term P^ et (M,r, R p ,p) in Eq. [TJ deserves special attention. In their original formal- 
ism, BG followed the common wisdom and assumed that the detection probability is mainly 
a function of the transit signal-to-noise ratio (SNR), which is usually defined by: 



SNR 



A 

cr 



(10) 



-6 - 



where n tT is the number of observed transits, A is the transit depth, and a is the photo- 
metric error. In BG's treatment, a, and therefore the SNR, strongly depend on the stellar 
magnitude, which is obviously the case in most surveys. 

However, unlike in most surveys, an important feature in the design of Gaia is a 'gating' 
mechanism, that will cause most bright stars, certainl y down to the 16 th magnitude, to be 
measured with a precision of 0.001 mag, more or less (IJordi et al.ll20fOI ). Thus, the SNR in 
our case is mainly a matter of the number of observations that occurred during transits. 

BG assumed that a planetary transit can be detected whenever the transit SNR exceeds 
a threshold. However, as opposed to high-cadence surveys, due to the small number of 
measurements that will sample the transits, Eq. [10] does not represent the complete detection 
problem in our case. Even if the SNR does exceed a threshold value, we still cannot regard 
the signal we found as a periodic one. In fact, in our assumed parameters, using the SNR 
naively, even a single observation in transit can be considered a transit detection with a SNR 
of 10. Obviously, we have to use another indicator, which mainly depends on the number 
of points observed in transit, and abandon the SNR figure of merit for our purposes (that 
would not be the case for smaller planets, though, where an even more elaborate indicator 
will probably be needed). To put it differently, we assume the problem is not limited by the 
error-bars of the individual measurements, but only by the scanning law and the temporal 
characteristics of the transit. 



von Braun et al.l ( 120091 ) calculated, for various ground surveys, the detection probability 
as a function of the period, which they dubbed the 'observational window function'. In our 
ideal case we assumed pure white noise, and no outliers. We further assumed a certain transit 
duration, and a minimum number of points in transit that would constitute a detection 
(depending on the detection approach used). We then calculated for each period, the fraction 
of configurations (namely, transit phases) which will result in detection, i.e., when the number 
of observations in transit will exceed the prescribed minimum. Once we had obtained the 
observational window function, we could integrate it over the required period range and 
obtain an estimate of the detection probability. 

The min imum number of observations in transit required for detection is not clear at 
the moment. iTingleyl ( 120111 ) introduced an algorithm, that requires a minimum number of 
seven to eight points in transit to secure detection. In an upcoming paper (Dzigan & Zucker, 
in preparation), we show that we can detect transits w ith five points in transit, or use our 
'Directed Follow-Up' approach (IDzigan fc Zuckerll201ll ) even for three points in transit. We 
therefore repeated our calculations here for a minimum number of three, five and seven 
points in transits. 
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We divided the entire sky into rectangular patches, 15° x 15° (apart for the poles, 
obviously, where we simply used the remaining circular patch), and applied the proper 
scanning law for each patch, assuming the scanning law for the central point as representing 
the entire patch. Fig. [T] shows a sample observational window function for one of the patches, 
for three cases of minimum points in transit (three, five and seven), and for a transit duration 
of 2 hours. This specific window function represents an area that Gaia is expected to visit 
70 times. This i s the average expected number of measurements over the entire mission 



fide Bruiinell2012h . 



For comparison we also present the window function for an area with 130 measurements 
in Fig. |2j The comparison shows that the detection probability depends strongly on the 
number of observations that the telescope will perform, during the mission lifetime. For 
example, the probability to sample a minimum of three transits (for an orbital period of 
3 days) increases from less than 30% for 70 measurements, to more than 60% in case the 
telescope should observe the star 130 times. 



Observational window function for 70 Gaia measurements 



(a=12 h 24 m 00 s ;8=-15°) 




P (days) 



Fig. 1. — A sample observational window function for a sky direction that Gaia is expected 
to observe 70 times, the average expected number of measurements over the entire mission. 
The detection probability is calculated for transit duration of 2 hours, for a minimum of 
three, five and seven observations in transit. 

Since we neglected the dependence of the window function on the SNR and therefore on 
the stellar characteristics, it remains mainly a function of the period and the scanning law. 
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Observational window function for 130 Gaia measurements 
(a=00 h 00 m 00 s ; 5=-45°) 




P (days) 

Fig. 2. — A sample observational window function for an area that Gaia is expected to 
observe 130 times. 

We can therefore take it out of the integral sign in Eq. (TJ which we calculate separately for 
each patch. We also divided the apparent magnitude range {G = 6 — 16) into 1-mag bins 
and treated each bin separately. Only at the end of the integration we multiplied the result 
with the detection probability we obtained from the window function. 

3. Results 

Table [1] presents the resulting expected yield of transiting HJs and VHJs in Gaia pho- 
tometry, down to stellar magnitudes 14 and 16, for which we expect Gaia to observe on the 
order of 1.5 x 10 5 and 6.8 x 10 5 stars, respectively. Obviously, the required minimum number 
of observations in transit strongly affects the results, as well as the assumed transit duration. 
In Fig. [3] we present the transiting planets yield for a minimum of five sampled transits, for 
transit duration of 2 hours, divided into apparent magnitude bins, from M = 6 to M = 16. 
Our results show that the Gaia photometry is expected to yield on the order of hundreds or 
thousands of new planets, depending on the detection strategy. 

Given the large amount of stars that Gaia will measure, we must estimate the false- 
alarm rate. The probability that 1 mmag purely white noise will produce a single 'transit-like' 
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Expected HJs and VHJs detections in Gaia photometry 
for a minimum of 5 sampled transits 

500 1 ■ ■ . . . ■ . . ■ . ■ 1 



400 



300 



200 



100 




6 7 8 9 10 11 12 13 14 15 16 

G (mag) 

Fig. 3. — The expected yield broken into apparent magnitude bins. This yield was calculated 
for a minimum of five sampled transits, and a transit duration of 2 hours. 



outlier (with magnitude of 0.005 — 0.015 mag) is ~ 3.5 x 10 -7 , which amounts to ~ 2 x 10~ 15 
for three outliers out of the average 70 measurements. This is a worst case estimate, as 
requiring the outliers to have a periodic pattern, even reduces this probability. Thus, under 
our nominal assumptions, it is obvious that the false-alarm rate is negligible. 

Nevertheless, the analysis might be complicated by stellar red noise. If the red noise 
doesn't possess a periodic or quasi-periodic nature, then the lo w-cadence sampl i ng sim ply 
renders it "white" , effectively reducing the SNR. According to McQuillan et al.l (120121 ) we 
estimate that roughly half of the stars have microvariability larger than 2 mmag. This 
results in an increased false alarm rate. However, in our upcoming paper (Dzigan & Zucker, 
in preparation) we show that our prioritization process in the Directed Follow-Up approach 
effectiviely eliminates them. 

In any case, Gaia will provide the astrometric and spectroscopic data needed for further 
classification. These data will help to exclude false-positives, such as background eclipsing 
binaries (BGEB), and to distinguish between periodic variability of stellar source and plan- 
etary transits. Thus, we can conclude that we don't foresee a significant false-positive rate, 
as long as we focus our analysis on HJs. 



Our results seem to differ considerably from those of iRobichonl (120021 ). We suspect that 
the main cause for the discrepancy is the different scanning law we used, which implied 
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50 — 200 o bservation s per star, with a mean of 70, compared to 100 — 300 with a mean of 



130, which iRobichonl used. This reflects changes in the mission design during the years that 
elapsed since 2002. Fig. [T] and Fig. [2] hint at the significant effect this change had on the 
observational window function. 

It is important to stress again that the analysis we present did not consider smaller 
planets. Their detection will be more difficult, and moreover, the required follow-up, either 
photometric or spectroscopic, will be more complicated. This topic will require a much more 
elaborate and careful analysis. 

The analysis we present in this Letter is a rough estimate that is based on general 
assumptions. We obviously made some approximations on the way, but at the level of 
accuracy needed at this stage we feel they are justified. The most important conclusion is 
that it will be worth while to develop detection algorithms, that will be tailored for Gaia 
photometry, and that will be incorporated into the pipeline. This may reduce the minimiun 
number of points in transit required for detection, which will immensely affect the yield. In 
addition, establishing a follow-up network that will be able to respond to alerts will also 
have a crucial effect, again, through this redu ction in the number of required observations 



in transit (e.g. IWyrzykowski &: Hodgkinll201ll ) 



A significant feature one can notice while examining Table [T] and Fig. [3] is the very 
strong dependence of the yield on the limiting apparent magnitude. Extending the analysis 
to fainter magnitudes will require introduction of the SNR into the analysis. This effort 
will be useless unless high-precision radial velocities of such faint targets will be feasible. 
Extremely large telescopes such as the E-ELT may allow this kind of observations. The 
enormous increase in the number of detected planets with Gaia, if fainter stars are considered, 
may serve as a justification to build high-precision radial velocity spectrographs for those 
telescopes. 

Gaia will undoubtedly revolutionize astronomy in many aspects. Nevertheless, its con- 
tribution to the field of transiting exoplanets is usually expected to be marginal. The ex- 
pected transiting planet yield is the key factor to establish whether this field will benefit 
considerably from Gaia. Usually, transit surveys focus on dense fields, to maximize the 
chances to detect transits and effectively use their high cadence. Gaia, on the other hand, 
will be an all-sky, low-cadence survey. This kind of surveys are usually considered irrelevant 
for transit searches. The estimate we present here shows that Gaia will have a valuable and 
significant contribution also in this field, mainly due to its high photometric precision, and 
in spite of its low cadence. 
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Table 1. Gala expected detections down to a limiting magnitude 





Minimum Number of 








Points in Transit 


G = 14 


G= 16 






3 


230 


999 


w = 1 hr 


5 


42 


178 




7 


7 


30 




3 


596 


2605 


w = 2 hr 


5 


209 


902 




7 


73 


310 




3 


720 


3191 


w = 3 hr 


5 


364 


1577 




7 


156 


669 



Note. — The expected yield of HJs and VHJs from Gaia photometry, for three different 
transit durations, down to a limiting apparent magnitudes of G = 14, and G = 16. 



